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On uniformly generating Latin squares 



r^ I M. Aryapoor* and E. S. Mahmoodian^ 



Abstract 



By simulating an ergodic Markov chain wliose stationary distribution is 
uniform over tlie space of n x n Latin squares, Mark T. Jacobson and Peter 
Matthews 0], have discussed elegant methods by which they generate Latin 

^j ■ squares with a uniform distribution (approximately). The central issue is 

the construction of "moves" that connect the squares. Most of their lengthy 
paper is to prove that the associated graph is indeed connected. We give a 

tG • short proof of this fact by using the concepts of Latin bitrades. 



1 Introduction and preliminaries 



^ ' A Latin square L of order n is an n x n array with entries chosen from an n-set 

N, e.g. {!,..., n}, in such a way that each element of N occurs precisely once 
CN ■ in each row and column of the array. A partial Latin square P of order n is an 

nx n array with entries chosen from an n-set N, in such a way that each element 
of N occurs at most once in each row and at most once in each column of the 
array. Hence there are cells in the aiTay that may be empty, but the positions that 
are filled have been so as to conform with the Latin property of aiTay. For ease of 
exposition, a partial Latin square T may be represented as a set of ordered triples: 
{(i, j; Tij) I where element Tij occurs in (nonempty) cell {i,j) of the array}. 

, . Let T be a partial Latin square and L a Latin square such that T C L. Then T is 

r> I called a Latin trade, if there exists a partial Latin square T* such that T* n T = 

j^ ■ and {L\T) U T* is a Latin square. We call T* a disjoint mate of T and the pair 

T = {T,T*) is called a Latin bitrade. The volume of a Latin bitrade is the 
number of its nonempty cells. A Latin bitrade of volume 4 which is unique (up 
to isomorphism), is said to be an intercalate. A bitrade T ~ {T,T*) may be 
viewed as a set of positive triples T and negative triples T*. 
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Example 1 The bitrade I = [I, I*), where 

I ^ {{i, J; a), {i, j';b), {i' , j;b), {i' , j'; a)}, 

I* = {{-i, j;b), {i, j'; a), {i' ,j; a), {i' ,j' lb)}, 

is an intercalate. Note that we must have i =/= i', j ^ j' and a =/= b. Usually such 
an intercalate is shown as 





j 


f 


i 


■ a^ ■ 


■ ba ■ 


i' 


■ ba- 


■ % ■ 



where the elements of I* are written as subscripts in the same array as I. 

For a recent survey on Latin bitrades see |2| and also |!5|. 

In p] the approach for generating Latin squares is based on the fact that an n x n 
Latin square is equivalent to an n x n x n contingency (proper) table in which 
each line sum equals 1 . They relax the nonnegativity condition on the table's cells, 
allowing "improper" tables that have a single — 1-cell. A simple set of moves 
connects this expanded space of tables [the diameter of the associated graph is 
bounded by 2(n — l)-^] and suggests a Markov chain whose subchain of proper 
tables has the desired uniform stationary distribution. By grouping these moves 
appropriately, they derive a class of moves that stay within the space of proper 
Latin squares. 

An improper Latin square is an n x n array such that each cell has a single 
symbol, except for one improper cell (in the improper row and column) which 
has three (the improper symbol appears there with a -1 coefficient). Each symbol 
appears exactly once in each row and in each column, except in the improper row 
(and also in the improper column) where one of the symbols appears twice as 
"positive" and once as "negative". An improper Latin square may be viewed as a 
set of n^ + 1 positive triples and one negative triple. 

Example 2 The following array is an improper Latin square of order 4. 



L 



c 


b 




d 


a 


b 


d 




a 


c 


d 


a + c- 


-b 


b 


b 


a 


b 




c 


d 



Using the notation of Latin bitrades, we may show this improper Latin square by 



L 



c 
b 

d a 
a 



b da 
d a c 
+ % b b 



The notion of ±l-move is introduced in 14]. Using the notation of Latin bitrades, 
a ±l-move means adding some appropriate intercalate to a given proper or im- 
proper Latin square such that the result is a proper or improper Latin square. If the 
added intercalate is the intercalate I in Example [T] the corresponding ±l-move 
is called a {{i,j;a), (i', j'; 6))-move. 

Example 3 By applying the ((1, 2; a), (3, 4; b))-move to the improper Latin square 
L in Example^ we obtain the following Lain square. 



L' 



c a d b 
b d a c 
d c b a 
abed 



Let G — (V, E) be a graph whose vertices are associated to S, the set of all 
proper and improper Latin squares of order n, and two vertices L and L' are 
adjacent if there is a ±l-move transferring L to L'. In the next section we state 
the results which prove that G is connected. This approach is developed from a 
linear algebraic approach to the concept of Latin bitrades, which is detailed in the 
references ||6), Q and Q. 



2 Connectivity of graph G 



In this section we prove that the graph G (defined in the last section) is connected. 
First we need a few lemmas. The first lemma states that an improper Latin square 
can be transferred into a proper Latin square using ±l-moves with changes only 
in two rows. 

Lemma 1 Suppose that we have the following improper Latin square 

3 



«2 



a +bs 



Then there is a sequence of (at most ^^^^) zLl-moves involving only rows ii and 
«2 which transfers A to a proper Latin square. 

Proof. It is easy to see that we can find the following cyclic pattern lying in rows 
ii and 12 of A (possibly after permuting some columns of A) 





i 


h 


J2 


h 


■ jr-1 


jr 


ii 


a +bs 


t 


u 


V 


Z 


S 


12 


s 


b 


t 


u 


y 


Z 



where i, w, ..., z ^ {s, a, 6} or r — 1 (i.e. t — s). Note that there is a similar 
pattern corresponding to a which has no intersection with the above pattern except 
in the jth column. Therefore one of these patterns is at most of length ^3^, and 
we may assume that r < ^^-5-^. We proceed by induction on r. If r = 1, then the 
{(*ii J; s), {12, ji', fo))-move produces a proper Latin square. Let r > I. Then the 
({«!) j;s)) (i2,jr;fe))-inove decreases r. ■ 

In the next lemma we show that one can swap a cycle lying in two rows using 
±l-moves. In ||9| it is called a cycle switch. 



Lemma 2 Suppose we have the following cyclic pattern in a (proper) Latin square 





ji 


h 


J3 


■ jr-1 


jr 


ii 


s 


t 


U 


y 


Z 


J2 


t 


u 


V 


Z 


s 



Then there is a sequence of (exactly r — 1) H-moves acting only on the entries 
shown above which transfers A to 





Jl 


J2 


J3 


jr~l 


jr 


h 


t 


U 


V 


Z 


S 


«2 


s 


t 


u 


y 


Z 



Proof. If r = 2 (i.e. u — s), then the ((ii, j'i;i), (12, J2; s))-move does the job. 
So let r > 3. Then the ((ii,ji;i), (12, j'2; s))-move transfers A to 





ii 


h 


J3 


■ jr-1 


jr 


n 


t 


s 


u 


V 


Z 


«2 


s 


U+ts 


V 


t 


s 



Now by applying the method in the proof of Lemma[T] this improper Latin square 
can be transferred to the desired Latin square. ■ 

The last lemma is a crucial lemma. It tells us that we can switch two entries in a 
row of an improper Latin square using a sequence of controlled ±l-moves. 

Lemma 3 Suppose for given s and t e {1, 2, . . . , n} we have the following im- 
proper Latin square: 



A 



where Z2 may be equal to i^. Then there is a sequence of (at most 2{n ~ l)j ±1- 
moves transferring this square to an improper {or proper) Latin square A' of the 
following form: 



A' 



where i' — 12 or 23, and the only possibly different entries of A and A' are entries 
in: (ii,ji), (Ji,j2) and those in rows 12 and i^. 





Jl 




h 


H 


S 




t 


«2 


a + 


bs 




J3 






s 





h 


h 


h 


t 


s 


i' 


^+ft 





Proof. We distinguish two cases. 

Case 1: i2 — is, i.e. in column J2 the symbol s appears in the improper row. So 
A has the following form: 



A 





ii 


J2 


ii 


s 


t 


12 


a +bs 


s 



Then the {{ii,ji;t), (12, J2; s))-move transfers A to: 



A' 





Ji 


]2 


h 


t 


s 


12 


a+bt 


t 



and we are done. 
Case 2: 12^ iz- 

It is easy to see that we can find the following cyclic pattern lying in rows i^ and 
13 of A (possibly after permuting some columns of A) 



A 





3\ 


Cl 


C2 


C3 


Cr-1 


Cr 


]2 




i\ 


s 


c 










t 




«2 


a +bs 


s 


U 


V 


X 


y 






h 


d 


u 


V 


w 


y 


z 


S 





where {m, w, w, . . . , x, y} fl {a, 6} = 0, but z G {a, 6}. Without loss of generality 
we assume that z = b. Therefore A has the following cyclic pattern 



A 





Ji 


Cl 


C2 


C3 


Cr-1 


Cr 


h 




h 


s 


c 










t 




12 


a +bs 


s 


u 


V 


X 


y 






h 


d 


u 


V 


w 


y 


b 


S 





Then the ((i2, Ji; s), (ii, ci; 5))-move transfers A to: 



Al 





h 


Cl 


C2 


C3 


Cr-1 


Cr 


J2 




ii 


b 


C+Sf, 










t 




12 


a 


b 


u 


V 


X 


y 






h 


d 


u 


V 


w 


y 


b 


S 





\fu = b (i.e. r = 1) then the ((is, ji; 6), (ii,ci; s))-move transfers Ai to: 





h 


Cl 


J2 




ii 


s 


c 


t 




12 


a 


b 






h 


d +bs 


s 


S 





which reduces the problem to Case 1 . So we assume that u ^b. Now the symbol 
b appears once more as a positive entry in column ci and another row, say i/^: 



Ax 





Ji 


Cl 


C2 


C3 


Cr-1 


Cr 


h 




H 


b 


c+s^ 










t 




i2 


a 


b 


U 


V 


X 


y 






h 


d 


u 


V 


w 


y 


b 


s 




ii 




b 















where ii ^ {«i, i2, is}. By Lemma[Tl with a sequence of ±l-moves only on rows 
ii and ii, we can obtain a proper Latin square A2. Using Lemma|2] a sequence 
of ±1 -moves transfers A2 to 



A. 





jl 


Cl 


C2 


C3 


Cr-1 


Cr 


32 


«2 


a 


u 


V 


w 


y 


b 




is 


d 


b 


u 


V 


X 


y 


s 



Now we can undo the sequence of ±l-moves on rows ii and 14 in ^3, to obtain 
the corresponding rows in Ai. The resuhing Latin square has the following pat- 
tern 





Ji 


Cl 


C2 


C3 


Cr~\ 


Cr 


k 




h 


b 


C+S5 










t 




i2 


a 


u 


V 


w 


y 


b 






h 


d 


b 


u 


V 


X 


y 


s 





With the ((is, ji; 6), (ii, Ci, s))-move, A^ can be transfered to 







ii 


Cl 


C2 


C3 


C,._l 


Cr 


k 






h 


s 


c 










t 




5 


12 


a 


u 


V 


w 


y 


b 








is 


d +bs 


s 


u 


V 


X 


y 


s 





and finally Case 1 finishes the proof. Note that in row ii, except the positions of 
s and t, other positions are unchanged. ■ 

Now we can prove that the graph G is connected. 

Theorem 1 Let S be the set of all proper or improper Latin squares of order n. 
Given two Latin squares of order n, there exists a sequence of ±l-moves that 
transfers one square into the other without leaving S. An upper bound on the 
length of the shortest such sequence is 2{n — 1)^. 

Proof. Suppose that A and B are two proper or improper Latin squares. Without 
loss of generality we can assume that A and B are proper (see Lemma [T]i. To 
prove the theorem, we proceed by induction on the number of identical rows of 
A and B. Suppose that the first fc — 1 rows of A and B are equal. We show that 
we can apply a sequence of ±l-moves to A to obtain a Latin square with the first 
k rows identical to the first k rows of B. If the fcth rows of A and B are equal 
then we are done. So suppose that they are not equal. In this case we can find the 
following patterns in A and B (s ^ a) 



A: 





ji 


k 


k 


s 


t 


k 


u 


s 





Ji 


J2 


k 


a 


s 


h 


■ b 


c 



Since A and B are proper Latin squares and have the same first fc — 1 rows, we 
musthaveii > k. Now the ((fc,J2; s), (ii,ji;i))-move transfers A to 





31 


J2 


k 


t 


S 


H 


u+st 


t 



which fixes the position of s in row k in both squares. 

If a — t, then we can find another entry in column ji of A which is equal to t 

and is not in the first k rows of A. So, by applying Lemma[Tl we can transfer the 

above (possibly improper) Latin square into a proper Latin square (using at most 

^^^^, ±l-moves) without changing the first k rows of A. 

If a y^ t, then we can find the following patterns in B and A 



A 





ii 


J3 


k 


t 


r 


h 


U +s^ 




«2 




t 



B = 





Ji 


h 


k 


a 


t 


ii 


■ b 




»2 




d ■ 



Since the first fc — 1 rows of A and B are the same we must have Z2 > fc. Therefore 
applying Lemma [3] interchanges t and r in row k without any other changes in 
row k and the first fc— 1 rows. Applying this process (at most n— 1 times) produces 
a (proper or improper) Latin square A' whose first fc rows are identical to those 
of B. Using Lemma [T] (and the fact that B is proper), we can transfer A' into a 
proper Latin square with a sequence of (at most ^^y^) ±l-moves. This finishes 
the proof by induction. In order to transfer A to i? we need to change n — 1 rows 
of A and for each row we need at most 2{n — 1)^, ±l-moves. Therefore with at 
most 2(n — 1)^, ±l-moves we can transfer Ato B. ■ 



Remark 1 Making moves "properly" 

In /¥/, they introduce moves that stay within the space of (proper) Latin squares. 
Such moves are called proper moves. Using Theorem\l\and a simple argument, 
they show that the space of (proper) Latin squares is connected under these proper 
moves. So we just explain what a proper move is in our notation. There are 
two kinds of proper moves, namely "two-rowed proper moves" and "three-rowed 
proper moves ". In order to define them, we first define the corresponding Latin 
bitrades. A tWO-rowed Latin bitrade is defined to be a Latin bitrade of the fol- 
lowing form: 





h 


h 


• jr-l 


ir 


h 


% 


be ■ 


yz 


Za 


«2 


ba 


% • 


^y 


az 



A three-rowed Latin bitrade is a Latin bitrade T with the following properties: 

1. T has exactly three nonempty rows, 

2. T is the sum of two-rowed Latin bitrades Ti and T2 such that there is at 
least one cell which is nonempty in both Ti and T2. 

Finally, a two-rowed proper move (resp. tinree-rowed proper movej means 
adding a two-rowed Latin bitrade (resp. three-rowed Latin bitrade) to a given 
Latin square provided that the result is still a Latin square. 

Another set of proper moves to connect the space of all Latin squares which is 
similar to the ones found by Jacobson and Matthews, but certainly found inde- 
pendently, appears in Arthur O. Pittenger fWj. Actually Pittenger's moves, corre- 
spond to special kinds of two-rowed and three-rowed moves, discussed above. 

Remark 2 The Markov chain introduced in id]/ is not known to be rapidly mix- 
ing (and thus does not have proven efficiency). Mark T Jacobson and Peter 
Matthews [4j state that: "in order to use either of our Markov chains to gen- 
erate almost-uniformly distributed Latin squares, we must know how rapidly the 
chain converges to the (uniform) stationary distribution. Of our two chains, we 
suspect that the "improper" one mixes more rapidly, in terms of real simulation 
time: executing a proper move takes time comparable to that needed to execute 
an equivalent sequence of ±l-moves; substituting an equal number of random 
-ill-moves seems likely to mix things up more." 

Remark 3 Randomly generating combinatorial objects is an important problem 
in combinatorics. It seems plausible to apply the ideas in this paper to attack the 



10 



same problem for some other combinatorial objects such as STS's. In fact one 
can define the notion of an improper STS see [1 ]. 
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